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1. Definitions Let P„ 
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l<i j'<n 



denote the n-by-n lower triangular matrices of binomial coefficients, Stirling partition 
numbers and Stirling cycle numbers respectively. For example. 
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We use V- and V for the falling and rising factorials respectively and we adopt the 
convention a binomial coefficient is zero if either of its parameters is negative. Let F^^r 
denote the 77,-by-n matrix ( (z — 1)^(*^-'^^^) j . Note that F^^r is a block matrix 

(c^.r o) with Gn,r lower triangular of size n — r. Let Hn^r denote the {n — r)-by-(n — r) 
lower triangular banded matrix ( (— 1)^*^-'^ (*Ci)'"^) • Note the entry of Hn^r 

V ^ ) l<i,j<n—r 

is ii i — j > r and can also be expressed as (^ll) (•~^)* Let Dn^r denote the diagonal 
matrix with the same diagonal 

/ 2 \ 

r - ^ H - ^ D 

^6,2 — ! -"6,2 — ' J^i 



12 6 
72 48 12 
V 480 360 120 20 / 



, \'i^)i<i<n-r, as Gn,r- For example. 
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2. Matrix Identities With the preceding definitions, we have the following four 
identities: 

/ 1 \ 

1 1 
2 1 



(1) 



n — 1 n / 



U. 



l<i,j<n 



l<i,j<n 



-1 



l<i,j<n 



diag((l,2,...,n)) 
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3. Canonical Forms Identity (1) puts P„ in bidiagonal form, which is good enough 
to obtain the Jordan form: say a matrix is near- Jordan if it has the form 

/A \ 

ci A 

C2 A 



Jn{^', Cl, . . . , Cn-l) 



\ 



Cn-l A / 



with all the q's nonzero. (A Jordan matrix is one with each q = 1.) Now J„(A; ci, , 
is similar to a Jordan matrix via 



, Cn-l) 



D J„(A; Cl,..., Cn-l) D ^ = Jn{\; 1, • • • , 1) 

where D = diag^(l, ci, CiC2, . . . , CiC2 ■ ■ ■ c„_i) j . Hence the Jordan form of P„ consists of 
a single block with A = 1. We can deduce the Jordan form of P„ mod p: reducing (1) 
mod p, the right side becomes block diagonal with as many p-hj-p near- Jordan blocks as 
will fit in an n-by-n matrix and (possibly) one of smaller size. Hence we have 



Theorem 1 The Jordan form of Pn mod p has precisely f^] blocks {all corresponding to 
the eigenvalue 1) and |_^J of them are of size p. 
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In particular, for n > p the minimal polynomial of Pn mod p is {x — ly 
Identity (2) shows 



Theorem 2 The matrix ( P.] ) = ( (— 1)* ■^l*} ) is a matrix of eigenvectors 

V / l<i,j<n \ } l<i,j<n 

for ((;)) 



Recall that two integer matrices A, B are equivalent if there exist unimodular matrices 
[/, V such that UAV = B. From (4), we have ( ° \'-) o) (^S''' J) = (^S" o) ^^^^^ 
the middle matrix in the product is -Fn.r- Combined with (3), we have 



Theorem 3 J„)^ is equivalent to the diagonal matrix diag{V, 2^, . . . , {n—rY, 0, . . . , 0). 



This is good enough Theorem 11.13, p. 30] to obtain the elementary divisors and hence 
the Smith normal form of (P„ — InY ■ 



4. Proofs After left multiplying by Sn, equating the (n + l,m + 1) entries and 
using the defining recurrence {^^^j^} = {^} + (m + j^}, identity (1) reduces to 

(m+i) ~ X]fc(fc){m} S 6.15, p. 265]. Similarly, identity (2) is equivalent to a simple 
variant of = [I] O B 6.16, p. 265] 

Equating the {n, m) entries, (3) is equivalent to 



^—^ h \ m I 1 r \ ^—^ 



k=m+r 



m j \ r I ^-^ \ r — 1 

k=m 



Both sides count the following set S of combinatorial objects: an object of is a parti- 
tion of [n] into cycles, m of which are colored black (say) and the rest are colored from 
{l,2,...,r} so that all r colors appear. The left side counts 5* by number of cycles: 
— partition [n] into k cycles; (^) — choose the black cycles; r!{^~™} — assign colors to 
the other k — m cycles. The right side counts S by total number of elements in the black 
cycles: — form a list of n — k elements from [n]; ("~^^^) — choose r — 1 of the n — k — 1 
"dividers" separating elements of the list, thereby obtaining a list of r nonempty lists each 
of which is canonically converted to cycles (for example, canonical might mean that cycles 
start at the left to right minima of a list), all cycles from the ith sublist getting color z; 
[m] — choose m black cycles on the remaining k elements. 
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Identity (4) is a consequence of the binomial theorem for rising factorials which implies 



that the inverse of = f f 

\ ■' / l<i.i<n—r V 



l<i,j<n—r V / l<ij<n—r 



5. Further Remarks The entry of (Pn — InY has a closed form, 
for which there is a simple combinatorial explanation. For a sequence c = (cj)j>o, set 
Pn{c) = (;-;)) . Then P„(c)P„(d) = P„(c*d) where c*d = (^^Lo ( •)c.t^n-.)„> 

If c counts labeled structures of some kind (c„ is the number of structures on [n]), then 
the r-fold convolution c* . . .*c counts r-sequences of structures whose labels partition [n\. 
In particular, for (the species of) nonempty sets c = (0, 1, 1,1,.. .) and c * . . . * c counts 
ordered partitions of [n] into r (unordered nonempty) blocks; the number of such ordered 
partitions is Hence, with c := (0, 1, 1,1,.. .) and d := ('^'{*})j>Q j {Pn — InY ~ 

P^{cr = P„(c * . . . * c) = P„(d) = (r!{7}(;.:;))^^. as asserted. 



Say an integer matrix is equivalent to its diagonal if it is equivalent to the matrix 
obtained by zeroing out all offdiagonal entries. Suppose a sequence c (such as ({*})j>Q) 
begins with r O's so that P„(c) has the block form (q„°(c) o) with (5„(c) square of size 
n — r. It follows from identities (3) and (4) above that the following statement is true. 
For c = ({r})i>o ' Qn{'^) equivalent to its diagonal. It seems the same statement is 
true for c = ([^])i>o would be interesting to have a proof. 
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